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An experimental investigation of a low 
Reynolds number turbulent boundary layer 
subject to an adverse pressure gradient 

By J. H. WATMUFF 


Summary 

This report covers a 12 month period from the beginning of the project. The 
project concerns an experimental study of a very low Reynolds number turbulent 
boundary layer subject to an adverse pressure gradient. The work is being per- 
formed in the Boundary Layer Wind Tunnel in the Fluid Mechanics Laboratory 
at NASA Ames Research Center in collaboration with Dr. R.V. Westphal. The 
aim is to obtain highly accurate mean-flow and turbulence measurements under 
conditions that can be closely related to the numerical simulations of Philippe 
Spalart for the purposes of CFD validation. It is expected that the experimental 
results will also serve as a useful contribution to the research literature in their 
own right. 

Much of the Boundary Layer Wind Tunnel has been completely rebuilt with 
a new wider contraction and working section which will improve compatibility 
with the simulations. A unique sophisticated high-speed computer controlled 3D 
probe traversing mechanism has also been integrated into the new test section. 
Construction of the tunnel and traverse has consumed a large fraction of the 
period covered by this report and is discussed in some detail. The hardware is 
now complete, and measurements are in progress. The mean-flow data indicate 
that a suitably two-dimensional base flow has been established. 

Automation of the probe positioning and data acquisition have led to a de- 
creased running time for total pressure measurements. However, the most sig- 
nificant benefits are expected to occur when using hot-wire probes. Calibrations 
can be performed automatically and there is no need to handle fragile probes 
when moving between measuring stations. Techniques are being developed which 
require sampling of the signals from moving hot-wire probes on the basis of their 
position in the flow. Measurements can be made in high intensity turbulence 
by flying probes upstream at high speed so that the relative magnitude of the 
turbulent velocity fluctuations are reduced. In regions where the turbulence in- 
tensity is not too large, the probe can also be repetitively scanned across very 
dense spatial grids in other directions. With this new technique, a complete 
profile can be measured in about 1/3 the time and with a spatial density about 
50 times that obtainable using a stationary probe. 
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1. Background 

1.1 CFD Simulations 

The numerical method used by Spalart (1988) is spectral in space and second- 
order in time. An important feature of the method is that there is no turbulence 
modelling. However, this requires a high grid density which restricts the simula- 
tions to low Reynolds numbers i.e. Re < 1500. Periodic spanwise and streamwise 
conditions are used and a multiple-scale procedure is applied to approximate slow 
streamwise growth. The key assumptions are that the streamwise evolution of 
the flow is slow and that the straining of the turbulence by the mean-flow can 
be neglected. Spalart has suggested that the assumption of small streamwise 
growth will cause the method to breakdown for large adverse pressure gradients. 

1.2 Relationship between the experiment and simulation 

There are three major requirements in the relationship between the experiment 
and simulation. Firstly, the Reynolds number of the experiment must match 
that of the simulation precisely. Secondly, the experiment and simulation also 
must have closely matched initial conditions. The simulation begins at the 
first station with an equilibrium boundary layer. A mildly favorable pressure 
gradient can be used to very closely approximate a “self-preserving layer” in 
the experiment, i.e. by careful experimental design the boundary layer can be 
maintained at almost constant thickness over some streamwise distance before 
being subjected to the adverse pressure gradient. Following a suggestion by 
Inman and Bradshaw (1981) the length of this region could be increased to allow 
any upstream trip effects to decay before the region of interest. Finally, accurate 
experimental pressure coefficient (C p ) measurements with high spatial resolution 
are required as an input for the simulation. A suitable flow configuration for 
the computations would be one in which the boundary layer experienced a non- 
dimensional pressure gradient ^ « 2 at a maximum Re « 1500. 

1.3 Preliminary experimental measurements 

The combined requirements of low Reynolds number and accurate pressure 
measurements precludes the use of water as a flow medium or the use of very low 
air velocities (e.g. 0.5 m/s) for the experiment. R.V. Westphal used momentum 
integral boundary layer calculations as a design tool to show that the Boundary 
Layer Wind Tunnel in its old form was a suitable facility for the experiment. 
Preliminary mean-flow measurements verified that this facility could duplicate 
the Reynolds number and pressure gradient requirements of the simulations. 
A flat test surface was employed for the boundary layer development, and the 
pressure gradient was imposed by a contoured upper control wall. The inlet flow 
velocity was 8.2 m/s. An inlet free-stream velocity around 8 m/s was considered 
to be close to the lower limit for accurate C p measurements, i.e. a total head of 
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0.15 inches of water with an accuracy of ±1%. The mean-flow two-dimensionality 
appeared to be good since the streamwise momentum balance was found to be 
within 10%. 

In the region of adverse pressure gradient at an Re « 1500 the reduced profile 
data showed a fl x « 0.7 which is a little low. A stronger adverse pressure gradient 
would be required for the current experiment. Also the inlet unit Reynolds 
number appeared to be a little high. A high inlet unit Reynolds number is 
undesirable since it shortens the physical distance over which the experiment 
matches the capabilities of the simulation. 

2. New Tunnel Contraction, Test Section and Traverse 

2.1 Tunnel Contraction 

When formulating the current experiment, it was desirable to use a lower inlet 
free-stream velocity which would reduce the inlet unit Reynolds number and also 
provide for better spatial and temporal resolution of hot-wire measurements. 
However, two factors had to be considered before a lower inlet velocity could 
be adopted. Firstly, the accuracy of the C p measurements would need to be 
maintained. Secondly, an increase in the total free-stream unsteadiness (from 
about 0.25% to around 0.5%) was observed at low tunnel speeds. 

The first restraint was relieved by purchasing the most accurate and sensitive 
commercially available pressure transducer. This allows the inlet free-stream 
velocity to be lowered to around 6 m/s (i.e. a total head around 0.08 inches of 
water) while still maintaining an accuracy of ±1% for the C p measurements. The 
second problem appeared to be related to the tunnel fan which was operating 
near the lower extremity of its characteristics. It was evident that a further drop 
in the fan speed using the original form of contraction and test section would 
lead to even more unsteadiness. Instead the flow area of the test section inlet 
was increased by 50% from 80cm x 20cm to 100cm x 24cm, allowing a lower 
inlet velocity but at a higher fan speed. This required designing a new 5:1 2D 
contraction to replace the original 7.5:1 3D version. 

Theoretical methods were used for selecting a new Boundary Layer Wind Tun- 
nel contraction design from a number of alternatives. The pressure distribution 
on the side walls was estimated using a potential flow method. The numerical 
solution of the stream function inside the contraction and inside upstream and 
downstream extensions of constant area was used to calculate the velocities on 
the boundary and the pressure distribution was obtained from the Bernouilli 
equation. The evolution of the contraction side wall laminar boundary layer 
subject to this pressure distribution was calculated using a momentum integral 
method. The effects of streamline curvature and convergence were neglected. 
The calculations showed that the new contraction design is not likely to cause 
separation. 
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2.2 High-speed probe traverse integral with test-section 

The numerical simulations usually have very high spatial resolution (e.g. 
Spalart (1988) used up to 10 7 grid points) but only for a limited number of 
time steps. This is in complete contrast to the experimental situation where 
almost unlimited temporal information is available but with very limited spatial 
content. One way of reducing this discrepancy is for the experimentalist to in- 
crease the spatial content of measurements. This dictates that multiple probes 
and/or more sophisticated traversing equipment be used. 

In view of the above considerations, a high-speed 3D computer-controlled 
probe traversing system has been built. The range of motion is 2.5m in the 
streamwise X-direction, 0.1m in the Y-direction normal to the wall and 0.5m 
in the spanwise Z-direction. Minimum top speeds of around 3 m/s in the X- 
direction and 1.5 m/s in the Y- and Z-directions were desired. These speeds 
can be reached over a short enough distance with accelerations of 2 g’s. Linear 
stepping motors satisfy the requirements for the Y- and Z-axes for accurate (0.05 
mm) high speed (1.5 m/s) positioning. The probe and sting are carried directly 
by the Y-axis motor. The Y-axis motor rides on a steel platen which is attached 
to the Z-axis motor, linear bearing and platen assembly. The Z-axis assembly 
is supported above the test wall within the working section by a gantry which 
spans the lm distance between the side walls. 

The requirements of the gantry are low mass (to minimize inertial loading), 
high stiffness (to minimize deflections), and small projected area (to minimize 
aerodynamic interference). Therefore, the gantry was constructed in the form 
of a lm wide truss using sections cut from thin sheets of carbon-fiber composite 
which were bonded together with epoxy resin. By itself the gantry weighs about 
2Kg which is less than 10% of the total movable mass. The total projected 
area of everything (i.e. the gantry, bearing rails, motors, platens etc.) located 
within the tunnel is equal to a one inch thick object spanning the tunnel side 
walls so the blockage is not excessive. Loading the gantry with a lOKg weight 
in the center of its one meter span results in a deflection of only 0.05mm. The 
ends of the gantry are fixed to carriages which move on linear bearings in the 
X-direction. 

The gantry is propelled in the X-direction by a brushless linear d.c. motor 
which reacts with a long magnet track to provide a maximum propulsive force 
of 80lbf. The motor is operated as a closed loop servo system with feedback 
provided by a linear quadrature encoder with a resolution of 10 micron. The 
maximum speed of the gantry has been tested at 3 m/s. Current linear stepping 
motor technology cannot match this performance. The two linear bearing rails 
are fixed to adapter plates which are bolted to a 1.2 x 3.0 m (4 x 10 foot) optics 
table. One of these adapter plates also supports the magnet track and linear 
encoder. The optics table provides an extremely rigid mounting platform for the 
entire assembly. Another significant benefit of the optics table is that it provides 
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an extremely flat reference surface for alignment purposes. 

Measurements are made in the boundary layer forming on a test plate which 
consists of a 0.375 ± 0.005 inch thick ground aluminum toolplate approximately 
lm wide by 2.1m in length. The plate is supported above the optics table by 
machined standoffs. A dial gauge attached to the traverse indicated that the 
deviations of the test-plate over the entire traversable X-Z plane were within 
±0.15 mm. The X-axis motor and twin-rail linear bearings are located beneath 
the test-plate, and rubber strips are attached to the side walls and to the plate 
to seal the slots that provide access for mounting the gantry. The side walls and 
flexible ceiling which complete the test-section are attached to a frame that is 
bolted to the optics table. Full access to the test-plate and traverse assembly 
can be obtained in minutes by releasing the bolts and hoisting the frame away 
with a crane. 

3. Benefits of the Traverse-Some New Techniques 

Automation of the high-speed 3D traverse reduces the time taken for experi- 
ments since repetitive control commands do not have to be performed manually. 
Although this leads to decreased running time for total pressure measurements, 
the most significant benefits occur when using hot-wire probes. 

Hot-wire probes can only be used with an acceptable accuracy when the mean 
velocity component is large compared to the velocity fluctuations. For a normal 
wire, the assumption of sensitivity only to the streamwise velocity fluctuations 
will become poorer as the fluctuations of the o'ther component normal to the 
wire increase (by simple vector addition). Crossed-wire probes have a theoret- 
ical velocity vector wedge-angle limit defined by the normals to each wire. In 
practice, the allowable flow angle is much less than this owing to wakes which are 
shed from the prongs. The boundary layer turbulence encountered in the region 
of adverse pressure gradient is expected to increase to such an extent that these 
sources of error will become significant. One way around this problem is to fly 
the probe in the upstream direction at high speed while sampling on the basis of 
the probe position. The superimposed bias velocity will reduce the magnitude 
of the turbulent velocity fluctuations relative to the increased velocity seen by 
the probe so that accurate measurements are possible. 

There is another, very distinct, advantage of flying the wire, i.e. measurements 
can be obtained on a very dense spatial grid. The time required to obtain 
converged data by repetitively scanning a dense grid can be more than an order 
of magnitude less than the time required to get the same data by traversing the 
probe to each point on the grid in succession and taking data with a stationary 
probe. At the low speeds of this experiment the sampling period must be as 
long as 100 seconds for adequate data convergence. In this case, the total time 
to obtain 30 points in a profile in the usual way with a stationary probe is 
of about 50 minutes duration. In regions where the turbulence intensity is 
not too high, the probe can be scanned up and down through the layer at 
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speeds that are quite high but not large enough relative to the mean flow to 
exceed the thresholds mentioned above. For example, with constant magnitude 
accelerations of alternating sign and a maximum speed limit of around 0.6 m/s, 
the oscillation frequency for traversing a distance of 32mm up and back down 
through a layer is about 4.5Hz. If samples are taken on the basis of the probe 
position and if 5,000 samples at each point are adequate for data convergence, 
then only 18 minutes is required to measure the entire profile compared to 50 
minutes in the case of the stationary probe. However, a more significant benefit 
is that there is no time penalty for increasing the spatial density of the grid. For 
example, if samples are taken for each 0.025mm increment of the probe position, 
then the spatial density of the grid can be increased by a factor of more than 40 
compared to the conventional method with a stationary probe. Also, spanwise 
derivatives and streamwise derivatives can be measured. 

Calibration of crossed-wire probes requires a mechanism in addition to the 
probe traverse. For example, in dynamic calibration schemes the small per- 
turbation sensitivity of the probe is determined directly by oscillating it back 
and forth continuously in a steady and uniform free-stream using a mechanical 
shaker. An alternative calibration scheme consists of imposing a number of flow 
angles on the probe by rotating it about an axis normal to the plane of the 
prongs. Intricate mechanisms are usually added to the traverse to perform this 
function. However, the same effect can be accomplished relative to the probe 
by moving it at high speed across a uniform free-stream. The advantage of a 
high speed 3D traverse is that either crossed-wire calibration scheme can be 
used without the need for additional hardware and without the need to handle 
probes. 

Hot-wire probes must be calibrated frequently in a uniform stream. Using a 
conventional single-axis traversing system in the adverse pressure gradient ex- 
periments would require the removal and reinstallation of fragile and expensive 
probes from the measurement region for the purpose of calibration. This would 
be time consuming and risk probe breakage. With a high-speed traversing sys- 
tem this could be performed almost instantly as well as minimizing the chances 
of probe damage. 

Since all the functions described above are under software control, the oppor- 
tunity exists to create an intelligent system. For example, hot-wire calibration 
drift is a particularly common and frustrating phenomenon during the course 
of an experiment. At regular intervals, the probe could be rapidly moved to 
a reference point in the free-stream to check for drift and another calibration 
performed (automatically) if some tolerance level was exceeded. 

4. Results to Date 

4.1 Incoming laminar boundary layer 

The boundary layer on the test plate has its origin upstream and, therefore, 
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it is influenced by the pressure gradients, wall curvature and lateral conver- 
gence within the contraction. Since the flow is accelerated by a contraction, 
the pressure gradients on the side walls are generally favourable. However, it 
can be shown that for a finite length contraction two regions of adverse pres- 
sure gradient exist on the walls in the vicinity of the entrance and the exit. 
Laminar boundary layers are especially sensitive to adverse pressure gradients 
and are prone to separation causing increased unsteadiness. Because of the low 
Reynolds number requirements of the current experiment, the quality of the 
incoming laminar boundary layer needs even more attention than usual. 

Laminar boundary layer profiles were measured with a total pressure probe 
at a station 8.3cm downstream of the contraction exit for a number of free- 
stream velocities. The pressure was constant along the test section for these 
measurements. All the velocity profiles were found to be very near the Blasius 
profile. In particular the boundary layer corresponding to the lowest velocity of 
U 0 = 6.3 m/s planned for the experiment has a shape factor H=^-=2.45 which 
is close to that of the Blasius profile (H=2.59). 

4.2 Selection of a transition device 

The low Reynolds number requirements of the experiment dictates that mea- 
surements be obtained close to the boundary layer tripping device. The effect 
of a simple cylindrical wire has been the subject of many studies, see Schlicht- 
ing (1979) p. 537 for a review. More recently there have been observations of 
spanwise irregularities in the boundary layers behind trip wires. Many workers 
have suggested that distributed three-dimensional roughness elements may be 
superior for transition purposes. However, the author is not aware of a system- 
atic parametric study that offers a reproducible alternative to a trip wire. The 
observed irregularities have varying strengths depending on the facility. There is 
evidence to suggest (see Bradshaw 1965) that the observations are more closely 
associated with wind tunnel screens rather than with trip wires. 

For the reasons outlined above, a simple wire was selected for boundary layer 
transition. The diameter and streamwise position for a number of wires were 
determined using the guidelines reported in Schlichting in conjunction with the 
laminar velocity profile measurements discussed earlier. Satisfactory results were 
obtained with 1.7mm and 2.0mm (diameter) wires when positioned near the 
contraction exit. However, when moved to a location 20cm downstream of the 
contraction, the 1.7mm wire failed to produce transition at all. Transition was 
intermittent with the 2.0mm wire at this location. These observations can be 
explained in terms of the adverse pressure gradient (not measured) near the 
exit of the contraction. Effective transition was obtained with a larger 2.4mm 
diameter wire at X=15cm. Measurements of Cj, momentum thickness and 
shape factor H indicated that a regular turbulent boundary layer is established 
behind the 2.4mm wire by X « 25cm corresponding to a momentum thickness 
Re « 450. 
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4.3 Mean velocity profiles with the pressure gradient 

The variation of C p with streamwise distance is shown in figure 1(a). The 
streamwise spacing (5cm) of the measurements appears to be adequate since 
estimates of dC p /dX obtained with four different numerical schemes are much 
the same as shown in figure 1(b). Mean velocity profiles were measured on the 
centerline (Z=0) with a 0.042 inch outside diameter total pressure probe which 
also functioned as a Preston tube for estimating the wall shear stress. The same 
probe was also used for the spanwise wall shear stress surveys. Mean pressures 
were determined over a 90 second period requiring about 45 minutes per profile 
and the data was acquired over a three day period. The probe was traversed 
under-computer control and once set running in the morning, the experiment 
proceeded without manual intervention throughout the day. 

A number of interesting features can be determined from the data set. A list 
of quantities derived from the profile data is given in table 1 and some of these 
have been plotted in figures 2(a)-(f). For 35cm < X < 55cm the layer thicknesses 
are approximately constant and the shape of the velocity profiles is similar as 
indicated by the shape factor. This suggests that the layer in this region is close 
to equilibrium. The profiles are shown in both wall- and outer-coordinates in 
both the region of the favorable and the adverse pressure gradients in figures 
3(a)-(d). Estimates of the skin friction coefficient Cf 0 obtained from the Preston 
tube have been used to plot the profiles shown in wall-coordinates. The pressure 
gradients are well within the allowable limits suggested by Patel (1965). The 
deviations from the log law indicate that uncertainties of up to 10% in wall shear 
stress may have to be accepted. 

4.4 Checks for two-dimensionality 

The results of two spanwise C/ surveys are shown in figure 4(a). The location 
of the upstream survey corresponds to the region where the layer is in equilib- 
rium. The survey spans approximately 40 boundary layer thicknesses and the 
variations are within ± 5% of the mean value. The downstream location in the 
region of adverse pressure gradient corresponds to Re « 1400 which is within the 
capabilities of Spalart’s simulations. The spanwise extent of the survey is about 
20 layer thicknesses, and again the variations are within ± 5% of the mean value. 
The momentum balance shown in figure 4(b) is within 10%. The spanwise Cj 
measurements and the momentum balance indicate that the flow is acceptably 
two-dimensional. 

5. Conclusions 

The boundary layer in the remodelled Boundary Layer Wind Tunnel is closely 
related to the capabilities of Spalart’s simulations. In particular, the layer ap- 
pears to be close to equilibrium just before application of the adverse pressure 
gradient which is a suitable initial condition for a simulation. The pressure gra- 
dient parameter /3 X « 2 at R e « 1650 which is about the upper limit of the 



Turbulent boundary layer subject to an adverse pressure gradient 


161 


numerical simulation capabilities. At this location, the distance downstream of 
the contraction is equal to the width of the test section. The spanwise variation 
of Cf is within ± 5% and the momentum balance is good. 

The high-speed computer-controlled 3D traverse has unique capabilities that 
will allow hot-wire measurements with a high spatial density to be obtained very 
quickly. Accurate measurements can be obtained on the “fly” in regions of high 
turbulence intensity that would be otherwise be subject to gross uncertainties 
using stationary probes. 
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TAB LEI. Streamwise Development of Mean Velocity Profiles with Pressure Gradient 


NOTE: 

1. Reference inlet velocity U # %6.$ m/s at X=Qcm. 

2. Boundary layer tr : p wire k at X=15cxn. 
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Figure 1(a) Variation of Pressure Coefficient C ? with streamwise distance, 
(b) Estimates of dC p /dX using four numerical schemes. 







X. cm x, cm 

W (f) 


Figure 2. Variation of quantites derived from profile data with streamwise distance. 

(a) Skin Friction Coefficient C/ (b) Displacement 6* and momentum 0 thicknesses 
(c) Reynolds number based on momentum thickness R $ (d) Boundary layer thickness 6 
(e) Shape Factor H=y. (f) Pressure gradient parameter 








Turbulent boundary layer subject to an adverse pressure gradient 



u 


Figure 3. Mean velocity profiles. 

(a) Region of favorable pressure gradient and (L) in wall coordinates, 
(c) Region of adverse pressure gradient and (d) in wall coordinates. 





Figure 4(a) Spanwise variation of C / at X=50cm and X=90cm 
(b) Momentum balance. 







